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1^ Summary. We described here our perception of complex systems, of how we feel 

I the different layers of description are an important part of a correct complex system 

simulation. We described a rough models categorization between rules based and 
i i law based, of how these categories handled the levels of descriptions or scales. We 

then described our fluid flow simulation, which combines different fineness of grain 
in a mixed approach of these categories. This simulation is built keeping in mind an 
ulterior use inside a more general aquatic ecosystem. 
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>^ 1 Introduction 

, . Not everyone agrees on what a complex system is, but many authors agree on 

some properties of these objects, among them the difficulty of their simulation. 
t-H We will describe here some of the causes we perceive of this difficulty, notably 

the intricate levels of description necessary to tackle their computer simulation 
. We will then explain how this notion guided our simulation, a fluid flow 
computation with solid obstacles, meant as a step toward a broader simulation 



CO 

of an aquatic ecosystem 

<N 

I>- 2 Simulation of complex systems 

2.1 Complexity vs reductionism 

•i-H 

As explained in [1], to study how mass works in a material system, dividing 
this system into smaller parts is a good method. Indeed, each of his subsystems 
is massive, and therefore the study, the reductionism, can continue. 

This is not so with living systems. If you divide a living system into smaller 
parts, the odds are good that all you reap is a heap of dead things. That's 
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because the life question of a live system is complex. This means that what 
is important is not so much the parts of the systems, nor the parts of these 
parts, but the functioning relations that exist between them. 

This is one of the two main reasons why one may want to integrate the 
multiple possible scales of description into a simulation. When you enquire 
about a complex question in a system, you need to choose carefully the needed 
levels of description, as you can't simplify them. Furthermore, these needed 
levels may change during the simulation, and it would be a fine thing if the 
simulation could adapt to these variations. 

Thus changing the scales of description during the simulation could be useful 
for the accuracy of the answers to complex questions regarding the system 
the simulation may provide . 

Then, there is the understanding of these answers. 

2.2 Clarity of the simulation 

Users of a simulation question it . Final users ponder about the future of 
the thing simulated in various circumstances, developers try to ascertain the 
validity of their model and of its implementation, but all use it with a purpose 
in mind. 

Choosing the right level of description is then important to give a useful 
answer. If the simulation is able to adapt its descriptions to what is needed by 
its user, lowering the noise and strengthening the signal, by choosing the right 
level (s) of description, it will be a better tool. For example in our application, 
a simulation of a fluid flow in an ecosystem, this help takes the form of hiding 
tiny perturbation and putting forward the main structures of the flow that 
emerged during the simulation. 

3 Methods for changing the scale in a simulation 
3.1 Law-based vs. rule-based models 

Classifying the various ways science can tackle problems is an arduous task. 
We will nonetheless distinguish two rough categories of models. 
Law-based models are the most used in science, most notably in physics. 
They are often continuous, especially in their handling of time and space, and 
based on a differential formulation whose resolution, ideally formal but often 
numerical, computes the values of state variables that describe the studied 
domain. Those methods are sometimes also called global or analytical. 

In rule-based models , the studied domain is discretized in a number of 
entities whose variations are computed with rules. There is therefore no longer 
a global description of the domain, nor is there a priori continuity. Cellular 
automata fall in this category of course , and so do objects/actors/agents 
. Those models have had a strong influence on game theory, and from there 
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directly on social models, and later on other domain through computer science 
for instance, at least by way of metaphors. Those models have other names 
depending on the domain where they are used, ranging from micro -analytical 
in sociology, to individual-based in life sciences or just simply local. 

Both kinds of models can be deterministic or stochastic. Finally, to blur 
the distinctions a bit more, models may include sub-parts falling in any of 
these categories. This is often the case with ecosystems for instance. 

3.2 Changing the scale in law-based models 

Accessing different levels of description in these models is often done through 
integration. Indeed, as said before, state function in these models are often 
continuous, and can therefore be integrated. New state functions are then val- 
ued or even built, on another domain and based on different phenomenological 
equations. For example, A. Bourgeat [4] describes fluid flows in porous milieus, 
where, from Navier-Stokes equations, through integration and the addition of 
an extra parameter, he builds a Darcy law. These changes of equations de- 
scription from one level to another alter sometimes drastically the linearity of 
the models and may lead to the introduction of new parameters that act as a 
memory of the local domain inside the global one. 

In a similar way to this example, the change of level of description in 
analytical models is often performed a priori, at the building of the model. 

3.3 Changing the scale in rule-based models 

Models based on rules offer a wider variety of ways of changing the levels of 
description. Indeed, local approaches are better designed to integrate partic- 
ularities of very different entities and their mutual influence, as is the case 
when entities of various scales interact. 

Cellular automata 

The first individual based computer science structures may have been cel- 
lular automata. If they were created by Stanislas Ulam, Von Neumann self- 
replicating automata may have been the foundation of their success [18]. Ulam 
himself already noticed that complex geometric shape could appear starting 
with only simple basic blocks. Von Neumann then Langton [10] expanded this 
work with self-replicating automata. 

If shapes and structures did appear in the course of these programs, it 
must be emphasized that it were users, and not the programs themselves, that 
perceive them. Crutchfleld [6] aimed at correcting that trend, by automating 
the detection of emergent structures. 

Detecting structures has therefore been tried, but reifying these structures, 
meaning automatically creating entities in the program that represent the 
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detected structures has not been tackled yet, as far as cellular automata are 
concerned. It could be that the constraint on its geometry and the inherent 
isotropy of the cellular automata are in this case a weakness. 

Ecology 

Since the beginning of the use of individual based models in ecology, the 
problem of handling the interactions between individuals and populations 
occurred [7] . The information transfers between individual was handled either 
statistically [5] or through the computing of action potential [14]. 

DAI uses 

Most software architectures designed to handle multiple levels of description 
are themselves hierarchical. They often have two levels, one fine grained and 
the other coarse grained. Communication between these two levels could be 
called decomposition and recomposition, as in [13]. 

In 1998, members of the RIVAGE project remarked that it was necessary 
in multi-agent simulations, to handle the emergent organizations, by associat- 
ing them with behaviors computed by the simulation [15]. Before that, were 
handled only border interactions between entities and groups [9]. 

This led in D. Servat PhD thesis to a hydrodynamic model incorporating 
in part these notions. In his Rivage application, water bowls individuals are 
able to aggregate in pools and rivulets. The individuals still exist in the bigger 
entities. The pros are that it enables their easily leaving the groups, the cons 
that it doesn't lighten in any way the burden of computing. Furthermore, 
these groups do not have any impact on the trajectories of the water bowls. 

4 Application to a fluid flow 
4.1 Ontological summary 

The fluid flows that constitute the ocean currents on the planet are the re- 
sult of an important number of vortexes of different scales. Turbulent move- 
ment can also be decomposed into vortexes, on scales going down to the near 
molecular. Viscosity then dissipates kinetic energy thus stopping the down- 
ward fractal aspect of these vortexes [12]. There are qualitatively important 
transfers of energy between these various scales of so different characteristic 
length. Representing these is a problem in classic modeling approaches. 

In classic, law based models, turbulent flows are described as a sum of a 
deterministic mean flow and of a fluctuating, probabilistic flow. These equa- 
tions (Navier-Stokes) are not linear, and space-time correlation terms must 
be introduced to compensate for that. These terms prevent any follow up of 
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the turbulent terms, and thus of the energy they transmit from one level to 
another. 

A pure law based approach is therefore not capable of a qualitative analysis 
of the transfer of energy between the different scales of a turbulent flow. A 
multi- level model, where multiple scales of vortexes would exist, and where 
they would be able to interact, would be a step in this qualitative direction. 

4.2 Treatment of multiple scales 

Fluid mechanic model and its structures 

There are a number of models used to describe fluid flows. The set we use here 
are based on a discretisation of the flow, and are called vortex methods [11]. 

In vortex methods, the flow is separated in a number of abstract particles, 
each being a local descriptor of the flow. These particles indicate the speed, 
vorticity etc ... of the flow where they are located. 

These particles are not fixed: they are conveyed by the fluid they describe. 
We find this model interesting as it is a local model, hence better able to deal 
with local heterogeneities. The values of the properties the particles describe 
are computed through the interactions between the particles, most notably 
through Biot-Savart formula. More details on this computation can be found 
in [3]. 

The vortex method we use is of 0(n 2 ) complexity. Finding ways of light- 
ening this calculus is therefore important. One lead is through making our 
model multi-scale, and only computing entities at the scale we need them. 
This is our second motivation for our using different levels of description. In 
order to have different levels of description, we will have to use an adapted 
description of the simulation entities. 

These entities come and go during the simulation, and thus we need a method 
to change the level of their description during the simulation, and not before- 
hand the way it is usually done. 

In our fluid flow, the main entities as we explained are vortexes. Not only 
do we therefore need to detect emerging vortexes by monitoring lower level 
vortexes particles, but also, as these vortexes aggregate among themselves to 
form even bigger vortexes, make this detection process iterative. Detecting the 
structures is not enough: we also need to create them in the simulation once 
they are detected. We must make these new entities live in the simulation, in- 
teracting with its various inhabitants (most notably particles, vortexes). They 
must evolve, whether it is growing or decaying to its possible disintegration. 

Let us now describe our recursive detection-creation/evolution-destruction 
cycle. 
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Detecting emergent vortexes among the vortex particles 

Structures are detected as clusters of particles sharing some properties. For 
vortexes these properties are spatial coordinates and rotation sense. As de- 
scribed in the following figure (figure 1), the process is: 

1. Delaunay triangulation of the particles 

2. Computation of a minimal spanning tree of this triangulation 

3. Edges that are too much longer than the average length of edges leading 
to the particles are removed. So are edges linking particles of opposite 
rotational. 

4. The convex hull of the remaining trees is computed 

5. An ellipse approximates the hull through a least square method. Geomet- 
ric constraints on the proportions of the ellipse must be satisfied or the 
vortex is not created. 

Further details on this process can be found in [16]. 



o Positive rotation 
n Negative rotation] 
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Minimal Spanning Tree 



Edges removing 



Fig. 1. Basic detection scheme 



Scale transfer : making simulation entities of the detected 
structures 

Detected structures are created in the simulation where they take the place 
of the particles whose interactions gave them birth. 
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The vortex structures are implemented through multiplicity automata [2]. 
These automata handle both the relations between higher level vortexes and 
the relations between them and the basic particles. 

Part of the relations between vortexes and their environment are handled 
through a method based on the eco-resolution model [8] , in which entities are 
described through a perception and combat metaphor. 




Fig. 2. Transducer managing the stability of a vortex 



The associated perceptions and actions are: 

1. Perceiving an intruder means being on a collision course with another vor- 
tex. Figure 3 sums up the various possibilities of interception by vortexes 
of opposed rotation and how each is translated in a perception. 

2. Attacking another vortex means sending it a message. 

3. Being attacked means receiving such a message. 

4. Fleeing means being destabilized: the vortex structure shrinks and creates 
particles on its border. (Figure 4). Too much flight can lead to the death 
of the structure, which is then decomposed in its basic particles. 

5. Getting satisfaction can mean two things. One is aggregating surround- 
ing particles of compatible vorticity. This calculation is done through a 
method close to the initial structure detection: Delaunay triangulation, 
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Fig. 3. Perception of an intruder 




Fig. 4. Destabilization 



spanning tree, removal of edges. Compacity criteria are then used to es- 
timate whether the tree should be added to the vortex and thus a new 
ellipse be computed or not. For instance in Figure 5, the particles on the 
lower left will be aggregated while those on top will not. The other is fus- 
ing with a nearby vortex of the same rotation sense, if the resultant vortex 
satisfies compacity geometric constraints (it mustn't be too stretched, as 
real fluid flow vortexes are not). 

The Eco-agent manages the stability of the structures. Its trajectory be- 
havior depends on adapted Biot-Savart formula. 
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Fig. 5. To aggregate or not to aggregate 



The described process is then iterated. New structures are detected and 
implemented, while others grow, shrink or disappear altogether. They move 
according to fluid mechanics laws. What remains to be seen is: how do they 
interact with solids? 

Further details on this process can be found in [16]. 

4.3 Interaction with solids 

Solids in a fluid flow simulation are a necessary evil, acting as obstacles and 
borders. Their presence requires a special treatment, but their absence would 
prevent the simulation of most interesting cases. They bring strong constraints 
in their interactions with the fluid, and as such introduce the energy that gives 
birth to the structures we're interested in. 

We manage structures using virtual particles. A solid close enough to real 
particles generates these particles symmetrically to its border (cf. Figure 6 
and Figure 7). The real particles perceive the virtual particles as particles in 
the Biot-Savart formula. The virtual particles on the other hand do not move. 
They are generated at each step if necessary, at the right place to repel the 
particles the obstacle perceived. 

Obstacles also generate virtual particles to try to repel structures. If it 
fails, obstacles aggress the vortex, thus making it shrink so as to deal with 
the interaction solid/fluid in a finer grain. 
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Fig. 7. Elliptic obstacle 
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Fig. 8. Simulation screenshot Red and blue : particles and vortexes Yellow disc : 
Obstacle In case of shades of grey, the disc is at the bottom left, surrounded by the 
fluid 



5 Conclusion 

We described here our perception of complex systems, of how we feel the 
different layers of description are an important part of a correct complex 
system simulation. We described a rough models categorization between rules 
based and law based, of how these categories handled the levels of descriptions 
or scales. We then described our fluid flow simulation, which combines different 
fineness of grain in a mixed approach of these categories. This simulation is 
built keeping in mind an ulterior use inside a more general aquatic ecosystem. 

Our result show so far an adequate real time handling of the interactions 
between fluids and solids. Other methods can statically simulate in better 
details this interaction, but it often requires a knowledge beforehand of the 
placement of the solids. We can deal with solids moving at random, or more 
interestingly computed by the simulation itself, and not known before the 
simulation begins. We hope it will enable us to simulate for example animals 
or dynamic obstacles, thus integrating our more global works on ecosystem 
simulation [17]. 
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